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Abstract 

The density of states for a particle moving in a random potential with a 
Gaussian correlator is calculated exactly using the functional integral tech- 
nique. It is achieved by expressing the functional degrees of freedom in terms 
of the spectral variables and the parameters of isospectral transformations of 
the potential. These transformations are given explicitly by the flows of the 
Korteweg-de Vries hierarchy which deform the potential leaving all its spectral 
properties invariant. Making use of conservation laws reduces the initial Feyn- 
man integral to a combination of quadratures which can be readily calculated. 

Different formulations of the problem are analyzed. 
PACS: 71.20.-b,71.20.Ad,71.25.Mg,71.55.Jv 
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The problem of a quantum particle moving in a random potential has been of interest for 
a long time. It has numerous applications to solid state physics particularly to the elec- 
tron structure of heavily doped semiconductors. Various modifications of the problem of the 
calculation of the Green's function and Density of States (DOS) averaged over the ensemble 
of random potentials arise (Thouless and Edwards f2|) in connection with the Anderson lo- 
calization || , in the physics of amorphous materials, particularly glasses |§] , in considerations 
of the conductivity in wires and mesoscopic objects ||, in the phenomena of broadening of 
resonance lines in solids |||7]], in a number of problems connected to the Quantum Hall effect 
H and many other problems related to disorder P-pl- 

In particular, the one-dimensional problem with a random potential having a white-noise 
correlator has been considered in a large number of works. The case when the potential is 
arbitrary but varies very slowly compared with the wave function of an electron was considered 
by Kane [|12|]. To calculate the DOS for the reverse situation with randomly distributed 
5-functional wells, Frish and Lloyd employed stochastic equations In the subsequent 

studies of the more general case of an arbitrary uncorrelated potential the methods of "minima 
counting" (Halperin and Lax [14]]) and "zeros counting" were used (Lifshits ||, Halperin |15|). 
Zittars and Langer formulated it in the Feynmann path integral language |16[] . That approach 
was further developed by Parisi who used the replica trick J4| and alternatively by Efetov by 
means of auxiliary Grassmannian fields ||17|| . Berezinski rigorously considered the problem 
by a diagram summation technique |]18[ . Luttinger and Waxier considered special potential 
configurations (5-functional scatterers) by the functional integral technique |L9|] . 

However, even for the ID problem with a white-noise correlator the average density of 
states i](E) has not been found exactly due to technical difficulties specific for the methods 
used. Most of the results obtained so far give only the "tail" behavior of the density of states, 
i.e., the DOS asymptotics at large negative energies E which is found to be ~ exp(— \E\ 3 ^ 2 ). 

A nonperturbative approach to the problem presented here uses essentially the ideas of 
Inverse Scattering Method (ISCM) for solving nonlinear evolution equations proposed in the 



celebrated works by Gardner, Greene, Kruskal and Miura p0|] , and Lax |2T| and developed in 
Refs. |22|-|2?f. Methodically, what is done in this work is a kind of inversion of this 'diamond 
of mathematical physics'. Here, the Schrodinger equation is to be solved and we use the 
hierarchy of Korteweg-de Vries (KdV) equations to find a parametrization of the random 



potential. 

DOS as a function of energy E averaged over the ensemble of random potentials with a 
Gaussian white noise correlator can be found from the functional integral (FI): 

(r](E)} = -\imIm[Tr(G E ^ 5 }} = 

71 8^0 

4=lim / D[U(x)]e- a3 I u ^ 2dx Im[TrG E - i5 {U(x)}], (1) 

71 Zi 8^0 J 

where Ge{U(x)} is the Green's function of a particle in the potential U(x), D[U(x)] is a 
functional integration measure, (...) means an ensemble average. The partition function Z = 
J D[U(x)]e~ a3 S u ^> 2dx is the normalization constant, the trace means the integral Tr[G] = 
Jf^, dxG(x, x' = x) and the constant a is the ensemble correlation length of the dimensionality 
[x] = {E' 1 / 2 ]. We use the units in which h 2 /m = 1. Due to the Gaussian weight only 
quadratically normalizable configurations enter the path integral ([[]): 

/+oo 
U 2 {x)dx < oo. (2) 
-oo 

Since the functions U(x) model the real impurity potentials, which are smooth, it is assumed 
that U(x) can be differentiated an appropriate amount of times. 

The Green's function G in (P corresponds to the eigenvalue problem for the Schrodinger 
operator L 

[E - L]G E (x,x') = S(x - x'), L^ k (x) = e k ^ k {x), 

t=-~ + U (x) (3) 

with the boundary conditions specified below. To average the DOS one has to solve these 
equations for every particular potential configuration U(x) entering the path integral and that 
appears to be the main difficulty in a straightforward approach. Let us note that the present 
problem ([]]) is similar to the one considered in the Random Matrix Theory (RMT) pq| . 

The boundary conditions for Eq.(|3|) are dictated by requirement @ and are typical for 
the scattering problem []. The particle is moving on a line in a square normalizable potential 
vanishing at infinity. Motion at E > is unbounded, and the spectrum remains continuous 



Actually, the stricter condition J|C7(x)|(l + |x[)dx < oo is required for rigorous validity of the 
following consideration. It means the potential vanishes rapidly at infinity. 



as for a free particle. Therefore the density of states at E > remains unaffected by such 
a potential and equals const x E~ l l 2 . Thus, we are interested in DOS for negative energies 
where all states are bounded and have discrete eigenvalues. So the boundary conditions for 
Eq.(|[) are assumed to be ip(±oo) = 0. 

Our method employs an accurate change of variables in the path integral ([|). The idea is 
to use the spectral data of a potential as a set of new variables instead of the potential itself. 
The spectral data alone cannot describe a potential configuration uniquely, because a lot of 
other configurations entering the integrand (|IJ) yield the same spectral properties. To account 
for all the potentials from the same isospectral family one needs to have an additional set of 
variables which describe all the spectrally equivalent configurations. The generators of the 
group of the isospectral transformations of the Schrodinger operator L are the Lax operators 
A 2 j + i obeying the commutation relations pi 



[d/dr 2j+1 + A 2j+1 , L] = 0, [A 2j+1 , A 2f+1 ] = 0, (4) 

where j is a positive integer and the "times" T 2 j+\ play the role of parameters and count the 
members of an isospectral family of potentials. The relations (f|) enforce the set nonlinear 
equations for U (x) (KdV hierarchy) 

U n = U x , U TS = -U XXX + 12UU X , ... (5) 

(the notation U a = dU/da is used). The role of parameters ti,t^,... is analogous to that of 
the angles of unitary transformations in the language of the RMT ||28|| . 

The essence of the ISCM theory (see, e.g., [22]- |[24|| ) which is important in the present 



context can be reformulated in the following rigorous statements: 

(i) Each square normalizable potential configuration, provided that condition (|]) is fulfilled 
sec footnote 1) can be characterized by an integer number of bound states N > 0, with the 

set of energies £j = k 2 , 1 < i < N, 

(ii) Each solution of the KdV hierarchy consists of a soliton component U s controlled by 
the finite number of discrete eigenvalues of the problem (J3|) and a radiation component U r 



2 This requirement is important for the establishment of analytical properties of solutions of (2) 
that used in the ISCM @. 



described by a continuous set of variables related to the reflection coefficient R. The soliton 
and radiation degrees of freedom are independent. 

(iii) The possibility for a configuration to evolve according to the KdV flows (|^) gives rise 
to a corresponding volume of isospectral transformations associated with it. The total weight 
of the configurations with definite set {ei, N} in the ensemble is proportional to this volume. 

Each of the potential configurations that enter our functional integral viewed as a solution 
of the KdV hierarchy fl5|) and it can be obtained from the trivial configuration U(x) = by 



means of the Backhand transformations or by Hirota's direct method [^5 . 

In view of the above classification (i)-(iii), the two following formulations of the problem 
of calculation of the average DOS are natural: 

A. The total number of bound states allowed for a given ensemble is fixed N{U(x)} = 
N = const for all particular realizations of the random potential. In that case, one should 
calculate the total level density i](E). 

B. The total number of discrete levels is not fixed by physical conditions, and calculation 
of the "reduced" DOS, i.e., the density per a length unit, rf{E) = dr)(E)/dx, is reasonable. 
In this case, one should keep the total "concentration" of bound states on the line c fixed. 

We shall consider both cases, implying that a particular choice should follow from the 
physical requirements for a given system under consideration. 

We evaluate the functional integral (1) using the exact parametrization of the potential 
viewed as a solution of the KdV hierarchy (5), that is given by the complete set of new 
variables [p3| , [26| ,p^[ 



U{x) = U(x, £j,{V(k)} , 6i,{Q{k)}) (6) 

actions angles 

where Sj < are the energies of the "solitons" (bound states), 9j are soliton phases re- 
lated to the soliton positions Xj as 0j = \2ej\ 1 ^ 2 Xj=J2m=o\^£j\ m+1 ^ 2 T2m+i, and the continuous 
set Q(k) defined at positive energies e = 2k 2 describes the isospectral transformations of 
the radiation part. The radiational actions set V(k) is related to the reflection coefficient: 
V(k)—— —log(l— \R(k)\ 2 ). As the KdV hierarchy is a fully integrable Hamiltonian system [23 



Eq.(^) is the canonical transformation from the potential variables U(x) with given number 
N of bound states to the action-angle type variables p6] . The canonical Hamiltonian struc- 



ture of the KdV hierarchy is governed by the standard Poisson brackets for the new variables 



{£j,0f} = Sjji, {P(k), Q(k')} = S(k — k') with all other brackets zero [f27|]. The Jacobean 
of the canonical transformation U(x) — > {ej,9j,V(k), Qj} is unity and the measure is 
just a product of differentials of the new variables. This transformation is correct because 
only smooth configurations are to be counted in the Feynmann integral and the problem with 
discontinuous paths typical in the field theories is irrelevant in the present context |K| . Thus 
we can express the functional measure in ([!]) in terms the new variables as follows 

j D[U(x)]<$>{U(x)} = [j[dU(x)${U(x)} = (7) 

J J X 

oo 

n ds 3 d9 3 D[V(k)]D[Q(k)MU(x, e jt 9 3 , {V, Q})}, 

3=1 

where $ is the integrand of Eq. QT|) . Let us stress that while we integrate over parameters of 
the evolution of the KdV equations we count all but the static potential configurations, and 
no real time evolution is involved in the problem at all. 

The next advantage of the soliton language is the possibility to rewrite the Gaussian weight 
in Eq. ([]]) compactly using the second conservation law for the KdV hierarchy |24| 
7 N 4 7 2k 2 

/ U{xfdx = T, J|2ei|i + / —logil - \R{k)\ 2 )dk. (8) 
J 3 J IT 

-oo * 

First we consider the case A when the total number iV of bound states is fixed. The 
expression flU) with the account for read as follows 

In 1 



x 



(V(E)} = = ^ J D[V{k)\D[Q{k)\ J Udside, 

p 2 fOO N 

exp[-4a 3 -a 3 / kdkVik))] Y S(E - e*), (9a) 

a 3 jo ~t 



with 



i 

4 poo 

xexp[ — a 3 y |2£j 3/2 -a 3 / kV(k)dk}. (9b) 
3 i Jo 

The (A^!) _1 factor appears due to symmetry reasons. The integration over the "radiation" 
variables {V(k), Qtk)} corresponds to the continuous spectrum of ISCM and is still continual. 
It is important that when expressed through the conserving quantities of the KdV equations, 
the integrands in (|9|) depend only on the action variables {Ei,V(k)}. 

The functional integration over V(k) and Q(k) factorizes giving the same functional con- 
stant / D[V(k)] D[Q{k)] exp[— a 3 / °° kV(k)dk] for both the numerator I N and the denominator 



Z N in the Eq.(^). It is the manifestation of the fact discussed above that the radiation com- 
ponent of the potential U r (x) does not yield bound states and does not affect the "soliton 
energies" £j. Hence it does not influence the DOS for E<0. To define the divergent inte- 
grations over the soliton phases Q{ we introduce the length of the system to play the role of 
a cutting parameter. In order to do so we rewrite the remaining integrals in the Eq.@ in 
terms of the soliton coordinates X* = 0j/\2eA 1 ^ 2 112 4)1. The Jacobean for this transformation is 



YliLi y2|£j| and we obtain from @ 

OO -,00 Jlft .art q 

(V(E)) = " 00 "°° - 00 (10) 

U?U dei\2ei\* J dX % ]exp[-^Y,tx\tei\* 12 ] 
-00 

It is interesting that similar expressions are normally obtained in the Gaussian ensemble of 
random matrices (except for the change e 2 — > e^ 2 ) |p8|] . However, the level repulsion is absent 
here contrary to the RMT, where it usually appears in the form of the factor Yli^j \&% — £j\ 
P%| . This property is specific for the the infinite length line, where one can have a potential 
configuration with two bound states with an arbitrary small difference between corresponding 
eigenvalues. Now all the integrals over the isospectral variables Xj cancel, and Eq. fllPf) is 
reduced to the combination of quadratures which gives the final result for the local DOS 



Case A : (rj(E)) = 4a 3 X ' ME\e-^ 2E ^ 2 (11) 



We turn now to the case B, where the number of bound states is not fixed from the very 
beginning. Due to the above classification every potential configuration belongs to one of the 
classes which differ by the number of the solitons N . And instead of (|S|) we can write the 
following expression for the averaged DOS 

= D I + D 1 I 1 + ... + D N I N + ... 
WV ;/ D Z + D 1 Z 1 + ... + D N Z N + ... V ; 

with the integrals In,Zn given by Eq.(^) 0. The dimensional functional constants Dm corre- 
spond to the relative probabilities for the N-soliton configurations to occur. To evaluate the 
expression (|1^) we recall that creation of a new N + l'th soliton in addition to N existing 
ones gives rise to an additional degree of freedom. It is related to "translation" of this extra 



Note the formal similarity between present classification and those arising in the nonperturbative 
field theories in discussions of the multi-instanton configurations [31]. 



potential well along the X axis. It enlarges the weight of iV+l-soliton potential configurations 
by a formally divergent factor / dX. Thus the configurations containing the number of soli- 
tons less than a given integer N have a zero measure in the ensemble compared to iV-soliton 
configurations, and Eq. fll^) reduces to 

{v{E)) = i im £^L = i im h (13) 

N^oo DnZn N^qo Zjy 



The fact that only the terms with maximum N survive in Eq. (|T2D can be proved explicitly by 
considering a discretized version of FI. The probabilities for iV-soliton potentials to appear 
are evaluated with the aid of the expression for the number of bound states (solitons) iV = 



\ Iu<o y\U(x)\dx valid for large N which we need in Eqs.(^J,|T3p. 

The total average DOS is related to the reduced DOS (i.e., per unit length) rj'(E) by the 
obvious relation i](E) = J dxrj'(E). The reduced DOS is independent of x due to translational 
invariance. After the simple transformations we obtain the result 



Case B : (rj'(E)) = Aa'cME^^ 2 ^ 2 (14) 



The total concentration of the bound states on the line, c, mentioned above, appears here as 
a natural limit N/ J dx — > c at very large N. 

Note that both the expressions (|TTD and (0) for the DOS correspond to the "local" density 
which is connected to the so called "cumulative" density of states J\f(E) by the relation 
r](E) = dM{E)/dE [0U]. The value N{E) is the number of levels with energy less than 
a given value E. Expressions (|TT]) and fll4l) are similar in form to those obtained elsewhere 
T3| . |T4|| for the asymptotic of the DOS in the large \E\ limit (band tails). However in the 



present approach no approximation has been used, and therefore the result is valid for all 
energies. 

To summarize, in this Letter we presented a nonperturbative method to calculate the DOS 
in a random potential. Using the exact classifications of the potential configurations in the 
functional integral in terms of the inverse scattering method, we obtain an analytical result for 
the DOS random potential with a white-noise correlator. No restriction by particular energy 
regions are needed within this approach. 

This method can be applied, though with some modifications, to a calculation of the 
averaged Green's function (G(x,x')) (two-point correlator). This will be presented in a forth- 
coming paper. 



Concluding, we note that the present approach is based on a deep connection between the 
ID Schrodinger equation and the exactly solvable nonlinear KdV hierarchy and hence can be 
directly applied to ID problems only But the method of using isospectral transformations 
to fulfill the averaging over the field ensemble looks quite fruitful and we believe that it can 
be done in a number of related problems in more than one dimensions in particular for the 
calculation of the DOS in a 2D random potential in the presence of magnetic field. It can also 
be used in random matrix theory. Work in these directions is now in progress. 

Comment added after submission: After our paper has been submitted to this 
database we were made known that a similar approach to the problem was discussed ear- 
lier by B.N.Shalaev in Fiz. Tverd. Tela (Leningrad) 32 3586 (December 1990), [Sov. Phys. 
Solid State 32 (12) 2079 (1990)] (Ref. [j3"3"] ). However, important aspects of considerations 
and the results of the two works are different. We are grateful to R.H. McKenzie for making 
us known about the work by B.N.Shalaev. 
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